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ABSTRACT 
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This paper proposes a novel radar system, namely step-frequency 
with compressive sampling (SFR-CS), that achieves high target range 
and speed resolution using significantly smaller bandwidth than tra- 
ditional step-frequency radar. This bandwidth reduction is accom- 
plished by employing compressive sampling ideas and exploiting the 
sparseness of targets in the range-speed space. 

Index Terms — step frequency radar, compressive sampling, nar- 
rowband radar 



1. INTRODUCTION 

Since the advent of radar systems much of the efforts have been de- 
voted to increasing radar range resolution. The relationship between 
range resolution and signal bandwidth is given by AR = where 
AR denotes range resolution, c is the speed of light and B is the 
bandwidth of the signal being used. Hence, wideband radar systems 
can achieve higher resolution than their narrow-band counterparts. 
However, wideband signals correspond to short pulses that experi- 
ence low signal-to-noise ratio (SNR) at the receiver. Further, they 
require high speed A/Ds, and fast processors 1 1 ]. Step-frequency 
radar (SFR) achieves high range resolution without sharing the dis- 
advantages of wideband systems. SFR transmits several narrowband 
pulses at different frequencies. The frequency remains constant dur- 
ing each pulse but increases in steps of Af between consecutive 
pulses. Thus, while its instantaneous bandwidth is narrow, the SFR 
system has a large effective bandwidth. Conventional step-frequency 
radars obtain one sample from each received pulse and then apply 
an inverse discrete Fourier transform (IDFT) on the phase detector 
output sequence for detection. Since the IDFT resolution increases 
with the number of transmitted pulses, SFR requires a large number 
of transmit pulses, or equivalently, large effective bandwidth. 

In this paper, we propose a step-frequency radar with compressed 
sampling, that assuming the existence of a small number of targets 
exploits the sparseness of targets in the range-speed space. The ap- 
plication of compressive sampling to narrow-band radar systems was 
recently investigated in |2]-|6|. A CS -based data acquisition and 
imaging method was proposed in 1 7 ] for stepped-frequency continuous- 
wave ground penetrating radars (SFCW-GPRs). In |8], CS-based 
step frequency was applied on through-the-wall radar imaging (TWRI). 
In 0,[ 8 ] the authors have assumed stationary targets and have shown 
that the the CS approach can provide a high-quality radar image us- 
ing much fewer data samples than conventional methods. Unlike 
1 7 1,| 8 ], the work in this paper explores joint target range and speed 
estimation based on compressive sampling, and proposes a radar 
system with reduced effective bandwidth as compared to traditional 
SFR systems. 



Compressive sampling (CS) |9|-|11] has received considerable at- 
tention recently, and has been applied successfully in diverse fields, 
e.g., image processing [12| and wireless communications |13|[14|. 
The theory of CS states that a K- sparse signal x of length M can 
be recovered exactly with high probability from G(K log M) mea- 
surements via £i -optimization. Let \l> denote the basis matrix that 
spans this sparse space, and let <E> denote a measurement matrix. The 
convex optimization problem arising in CS is formulated as follows: 
min ||s||i, s.t. y = <I>x = <I>\l/s = 0s where s is a sparse vec- 
tor with K principal elements and the remaining elements can be 
ignored; $ is an iV x M matrix with N <C M, that is incoherent 
with It has been shown that two properties govern the design of 
a stable measurement matrix: restricted isometry property and in- 
coherence property ifTTl . A if -sparse signal x of length-M can be 
recovered from N < M samples provided N > K and satisfies 

1 — e < ^IMlH 2 — ^ ~^ 6 wnere v * s an ar °it rar y if -sparse signal and 
e > 0. This property is referred to as Restricted Isometry Property 
(RIP). The incoherence property suggests that the rows {(/>;} of 
should be incoherent with the columns {^j} of * r - 

Let us consider an SFR system that transmits N pulses and waits 
for echoes of all pulses to return before it starts any processing. The 
frequency of k th pulse equals 



f k = f + kAf 



(1) 



where fo is the starting frequency and fce{0,l,2,...,iV— 1}. The 
k th transmit pulse is of the form rect(t) e ~ t27rfkt . 

The signal reflected by a target at distance R moving with speed 



2R 

r(t) — rect(t ) e 



-i27T/ fe (t-f (R+vkT)) 



(2) 



where c is the speed of light, T is the pulse repetition interval (PRI). 
Here we assume that v is small enough to be considered constant 
within the pulse interval. 

As we assume each pulse to be narrowband, we can ignore the 
delay of in the signal envelope, and just consider the phase shift 

(3) 



for extracting the information on targets. Therefore we have 
r(t)=rect(t) ^M^R+vkT)) 
and the output of phase detector is of the form 

y k = e ^Uo + kAf)^(R+vkT) 

The exponent of (|4} can be written as 



(4) 
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4tt , „ n 2AfR , n 2f vT 7 n 2kAfvT , _ 
7^ — f R + 2tt — k + 2tv— k + 2tt k (5) 

c c c c 

The first term in $5$ represents a constant phase shift due to the start- 
ing frequency, while the second term represents the phase shift due 



to frequency offset of the k— th pulse. The maximum unambigu- 
ous range and range resolution for step-frequency radar are given by 
R u = 2KJ anc * — 2jva/ res P ec tively. Here NAf is the total 
effective bandwidth of the signal over N pulses. Targets which are 
at distance R > R u will be seen by the system to be at distance 
R — R u . 

The third term of gives the Doppler frequency shift experi- 
enced by the signal due to the target speed v. The fourth term of 
(|5]) represents the frequency spread due to target speed. This has the 
effect of spreading the energy of the main lobe at the target position. 



3. THE PROPOSED APPROACH 

Let us take the transmitter and receiver to be co-located and em- 
ploy N pulses for estimating the range and speed of targets. Let us 
discretize the range space as [R±, . . . , Rm], and the speed space as 
[vi, . . . , vl]. The whole target scene can be described using M x L 
grid points in the range- speed plane. The range and speed spaces 



discretization steps are AR 



and Av 



L-l 



respec- 



tively. We assume that the targets can be present only on the grid 
points. By representing the target scene as a matrix S of size M x L, 
equation © becomes 

M L 

Vk=Y,Y,^ lk ^ Rm+vlkT) -S{rn,l)+w k (6) 

m=l 1=1 

where 

J a reflectivity of target present at (R 

I if target is absent at ((m — 1)L + l) th grid point 



(7) 



and Wk represents zero-mean white noise. 

Putting the outputs of phase of the phase detector, i.e., yk, k 
1, N in vector y, we get 

y — \l/s + w 



(8) 



where s = [si,S2, ...,sml] T , w represents white zero-mean mea- 
surement noise, and the elements of matrix \l/ equal 

i/>(k, (m - 1)L + Z) = e ^<fo^A f )i(R m+Vl kT) (9) 

for k = 1,...,JV, m = 1,...,M, I = 1,...,L. We can think of the 
basis matrix \l> as being a stack of column vectors {SE^}^ -1 , i.e., 



* = ( * | *1 |-..| *ML-1 



(10) 



where each is of size N xl containing the phase detector outputs 
for all N pulses corresponding to the phase shift due to a target lo- 
cated at the i th grid point. Thus, \l/ accounts for the phase shift of all 
possible combinations of range and speed. Taking the measurement 
matrix 4> to be an N x N identity matrix yields = 

Based on ([§]) we can recover s by applying the Dantzig selector 
to the convex problem ((T5)) 



s = min ||s||i s.t. ||0 H (r — 0s)||oo < 



(11) 



According to (T5), the sparse vector s can be recovered with very 
high probability if \i — (1 + t~ 1 )^/2 log Na 2 <Jmax, where t is a 
positive scalar, a m ax is the maximum norm of columns in the sens- 
ing matrix G and a 2 is the variance of the noise in ([6]). A lower 
bound is readily available, i.e., \i > log Na 2 a m ax- Also, \i 



should not be too large because in that case the trivial solution s = 
is obtained. Thus, we may set /i < HO^ r||oo. 

In conventional SFR systems, the IDFT algorithm requires the 
columns of the transform matrix to be orthogonal. The range res- 
olution in space depends on the frequency resolution in the Fourier 
domain. Therefore, these systems require N — M pulses in order 
to have a range resolution of ^ . For the proposed approach we can 
use N < M pulses and still achieve a range resolution of . 

For moving targets, the conventional IDFT method for estimat- 
ing range and speed observes a shift in the target positions (due to 
speed) and a spreading effect around the shifted position (due to the 
fourth term in equation ([5])). These effects degrade the receiver per- 
formance causing erroneous range estimation and sometimes miss- 
ing the target completely. In the proposed approach, since \l> has 
columns corresponding to all the possible range-speed combinations, 
the estimated results are comparatively more accurate. 

4. SIMULATION RESULTS 

Stationary targets - Our simulations use the following parameter 
values: fo — 1MHz, Af — lOKHz, number of grid points M — 
100. These values of Af give R u — 15 Km, and AR — 150 m. 
We assume that the stationary point targets are present on the grid 
points. 100 iterations of sparse target vectors are generated and the 
estimation accuracy is computed as the ratio of number of iterations 
for which the target ranges are correctly estimated to the total num- 
ber of iterations. 

The basis matrix \l/ of size N x M is generated according to 
equation ([9]). The measurement matrix <I> is an identity matrix I of 
size N x N. The optimization algorithm used to solve equation dTTb 
was obtained from fl6l . 

The number of pulses, N, controls the column correlation of 
the measurement matrix for a given value of M (number of grid 
points). Lowering the correlation between adjacent columns of ^ 
increases the isolation among the columns, which results in better 
range estimation. 

For the measurement matrix generated by using equation J9]), the 
adjacent column cross-correlation equals 



Rd) — 



1 1 



N i 



(12) 



Figure \T\ shows the effect of changing the column correlation 
on the estimation accuracy of the CS sensing matrix when the un- 
ambiguous range is divided into 100 grid points. Figure E] shows the 
accuracy of the CS detector in the presence of noise at different SNR 
values for the case in which only 5 targets are within the detectable 
range. The noise signal added at the received signal was Gaussian 
zero-mean with variance a%, where the variance a% changes with 
SNR. Figure[3]compares the performance of the CS detector with the 
conventional IDFT detector for N = 70 for a target scene contain- 
ing 5 stationary targets. As it can be seen, the CS detector performs 
better than the IDFT detector for all SNRs. Figures and [3] show 
that we can use N < M pulses to obtain AR = -j, provided that 
N — 0{K log N). This proves that we can use lower bandwidths 
in CS compared to conventional techniques of IDFT and accurately 
estimate the target parameters, when the targets are sparsely present 
in the range- speed space. 

Moving targets - The number of grid points in the range do- 
main is taken to be M — 40. L — 6 grid points were used for 
discretizing the speed axis. The carrier frequency is / = 10 8 Hz 
and the step frequency Af = 10 5 Hz. The ranges and speeds of 
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Fig. 1. Effect of changing column correlation on estimation accuracy Fig. 3. Accuracy comparison between CS and IDFT detectors for 
for 100 grid points for stationary target scenario N = 70 and 5 stationary targets. 
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Fig. 2. Accuracy of CS detector at different SNR values for 5 sta- 
tionary targets 



4 targets are generated randomly in each of 100 Monte Carlo runs. 
Out of 4 targets, two targets are placed on adjacent grid points on 
the range- speed space in each run. The accuracy of the detector 
is computed as the ratio of the number of runs for which all target 
ranges and speeds have been estimate accurately to the total number 
of runs. In Fig. 01 we show the detection accuracy of CS and IDFT 
methods for different values of SNR. For moving targets, the IDFT 
method requires speed compensation before performing IDFT. Since 
the target speed is unknown, we compensate the received signal with 
all possible speed and choose the one with the highest and sharpest 
IDFT output. As it can be seen, the proposed CS significantly im- 
proves the detection accuracy as compared to the IDFT method. The 
advantage of the CS approach is more obvious at low SNR. Figure 
[5] compares the detection accuracy of the CS and the IDFT methods 
for different number of pulses for SNR — 15dB. We can easily see 
that the proposed method requires much fewer pulses than the IDFT 
method to achieve the same accuracy level. For example, the CS 
approach requires 130 pulses to achieve detection accuracy of 0.95, 



while the IDFT method needs about 190 pulses in this particular case 
considered in our simulations. 
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Fig. 4. Accuracy comparison between CS and IDFT detectors for 
different values of SNR for N = 100 pulses (moving targets). 



One trade-off that is not apparent from the simulations is com- 
putation time. Convex optimization techniques have much higher 
computation cost as compared to the IDFT. The basis pursuit (BP) 
algorithm used in our simulations has a computation complexity of 
<3((ML) 3 ). Thus the processing speed of the receiver system may 
put a limit on the number of grid points ML that can be used and 
the range resolution AR. However there have been other algorithms 
like Orthogonal Matching Pursuit (OMP) which have computation 
complexities of O(NKML). A decoupled range-speed estimation 
approach along the lines of [ 17 1 could also be employed here to re- 
duce complexity. 
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Fig. 5. Accuracy comparison between CS and IDFT detectors for 
different number of pulses at SNR — 15dB (moving targets). 



5. CONCLUSION 

We have proposed a CS -based SFR system for joint range- speed 
estimation. It has been shown by our simulation results that the 
proposed CS approach can achieve high resolution while employ- 
ing lower effective bandwidth than traditional SFR systems. Unlike 
the IDFT method, the proposed approach does not suffer from range 
shift and range spreading around the shift positions caused by the 
movement of targets. 
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